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LINEAR  ESTIMATION  OF  HYPERSPECTRAL  MIXED  PIXEL 
COMPONENTS 

Edmundo  Simental  and  Elisa  C.  Gonzalez 
U.  S.  Army  Topographic  Engineering  Center 
Alexandria,  VA  22315-3864 


ABSTRACT 

This  paper  presents  a  method  to  determine  end  members  along  with 
their  relative  concentration  in  a  hyperspectral  mixed  pixel.  The 
method  is  modeled  as  a  linear  combination  of  the  end  members 
reflectance  spectra  and  a  library  of  spectral  prototypes.  This  method 
was  tested  with  431-band,  laboratory  controlled,  data  sets. 
Indications  are  that  this  method  can  be  extended  to  field  data 
observations  with  a  small  number  of  end  members  and  that  at  least  some 
error  source  can  be  identified  and  data  and  results  adjusted 
accordingly. 


1.  INTRODUCTION 


Linear  estimation  is  a  viable  method  for  analyzing  mixed  pixels 
under  the  simplified  assumption  of  accurate  reflectance  calibration, 
accurate  laboratory  spectral  prototypes,  and  dominant  linear  mixing. 

1 . 1  OBJECTIVE 

This  paper  addresses  a  method  for  estimating  the  concentration  of 
materials  in  a  mixed  pixel  using  laboratory  measured  hyperspectral 
data,  a  library  of  spectral  prototypes,  and  a  linear  model.  A  library 
of  spectral  prototypes  is  simply  a  collection  of  spectra  of  pure 
materials  or  objects  which  may  be  components  of  the  mixed  pixel.  The 
model  is  based  on  least  squares.  A  basic  requirement  is  that  the 
spectral  prototypes  be  linearly  independent.  The  library  can  contain 
several  hundred  prototypes  and  each  prototype  can  contain  several 
hundred  bands. 

1.2  HYPERSPECTRAL  MIXED  PIXEL  CONCEPT 

An  imaging  spectrometer  measures  percent  reflectance  or  radiometric 
values  of  a  pixel  at  various  frequencies.  A  plot  of  the  reflectance 
as  a  function  of  frequency  represents  the  spectral  response  of  the 
pixel.  A  pixel  is  said  to  be  mixed  if  its  measured  reflectance  is  a 
result  of  imaging  multiple  object  types  or  materials  within  the 
pixel's  footprint.  In  such  cases  the  measured  response  will  be  a 
combination  of  the  individual  responses  of  the  objects  or  materials. 
If  the  number  of  materials  in  the  pixel's  footprint  is  less  than  or 
equal  to  the  number  of  measurements,  and  if  the  combination  of  the 
responses  within  the  measuring  instrument  is  linear,  then  the  most 
probable  mix  of  components  may  be  estimated  by  a  linear  model  and  a 
library  of  spectral  prototypes  of  materials  suspected  in  the  mixed 
pixel . 


1.3  LINEAR  ESTIMATION  THEORY  CONCEPTS 


Estimation  is  defined  as  the  process  of  making  a  decision  or 
judgement  concerning  the  approximate  value  of  parameters  when  the 
decision  is  weighted  or  influenced  by  all  available  information  on 
those  parameters.  The  estimation  is  linear  when  the  estimation 
characteristics  are  additive  and  homogeneous . 1  In  this  particular  case 
additive  means  that  the  spectral  response  combination  is  the  sum  of 
the  product  of  the  individual  material  component  responses  and  a 
scalar  quantity  associated  with  the  particular  material.  Homogeneous 
means  that  the  scalar  quantity  is  dependent  on  the  material  and 
independent  of  the  frequency.  The  spectral  response  ’«*  at  frequency 
f0  is  composed  of  material  ’a'  and  material  '  b',  then  for  any  frequency 
*£ '  / 


Mf  =  k,( spectral  response  of  *  a'  at  frequency  ’ f  * )  + 

k2(  spectral  response  of  'b'  at  frequency  'f') 

where : 

«f  =  total  reflectance  at  frequency  'f' 

kj  =  scalar  constant  associated  with  material  'a' 

k2  =  scalar  constant  associated  with  material  'b' 

1.4  LINEAR  MODEL 

Assume  that  there  are  *p’  different  objects  in  a  spectrometer's 
footprint  and  that  'n'  discrete  frequencies  are  measured,  then  the 
output  reflectance  measured  at  any  frequency  'k'  is  given  by; 


w 


.-Jr 

‘  i=l 


k, 


where : 


G)k  =  total  output  reflectance  at  frequency  '  k'. 


rkfi  =  response  of  material  ' i'  at  frequency  'k' 
s^  =  coefficient  of  concentration  of  material  ' i' 

The  above  equation  in  matrix  form  is  given  by; 
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Ralph  Deutsch,  Estimation  Theory ,  Prentice  Hall 


The  above  expression  written  in  matrix  notation  is  given  by: 


r>=Rs 


(1) 


The  above  equation  represents  a  linear  statistical  measurement  model 
in  which  the  'R'  matrix,  (nxp)  is  .*  ’laboratory  prototype'  spectral 
library  of  reflectance  values  for  materials  that  may  occur  in  the 
pixel  footprint.  It  is  important  to  note  that  these  values  must  be 
obtained  separately  from  a  mixed  pixel  solution.  These  values  could 
be  obtained  from  measurements  in  the  laboratory  or  measurements  in  the 
field  under  the  same  conditions  as  those  employed  for  the  mixed  pixel. 
They  represent  the  spectrum  of  a  pure  material  or  objoct.  The  column 
matrix  'Q  '  (nxl)  represents  the  'n'  response  measured  by  the 
spectrometer.  The  column  matrix  's'  (pxl)  is  the  'coefficient  of 
concentration'  of  each  material  and  is  unknown.  The  problem  is  to 
estimate  this  parameter  from  the  'spectral  library'  and  from  the  'n' 
measurements.  Once  determined,  these  coefficients  can  be  converted  to 
represent  the  relative  quantity  of  each  material. 

1.5  LINEAR  MODEL  SOLUTION 

The  number  of  measurements  is  normally  expected  to  be  much  larger 
than  the  number  of  components  in  the  mixed  pixel  to  produce  an  over 
determined  system  of  linear  equations.  Methods  for  the  solution  of  an 
over  determined  system  of  linear  equations  are  many  and  are  well 
documented.  A  basic  least  square  solution  is  presented  here. 

Let, 

R*  =  transpose  matrix  of  R  in  equation  (1)  above  and  multiply  on  the 
left  both  sides  of  the  equation  by  it  to  obtain, 

Rc  Q  =  R*  Rs 

The  matrix  product  [R11  R]  produces  a  pxp  symmetric  matrix  and 
matrix  product  Rc  Q  produces  a  pxl  column  matrix. 

Let, 

m  =  inverse  of  [Re  R] 


(2) 


Multiply  on  the  left  both  sides  of  equation  (2)  by  m  to  obtain, 

m[Rc  Q  }  =  m[Rc  R]s  (3) 

By  definition,  the  matrix  product  m[R~  R]  is  the  identity  matrix  and 
equation  (3)  becomes 

m[Rc  Q]  =  s  (4) 

The  matrix  product  mfR1  Q]  is  pxl,  a  column  matrix,  and  is  the 
least  square  solution  vector  for  's'. 


2 .  LABORATORY  TEST 


The  development  of  an  understanding  of  potential  complications  to 
simplified  unmixing  models  and  their  ramifications  may  be  best 
understood  by  working  with  laboratory  controlled  data. 

2.1  EXPERIMENT  DESIGN 

A  laboratory  controlled  experiment  was  designed  to  assess  the 
validity  of  the  developed  model.  The  experiment  consisted  of 
obtaining  a  prototype  'spectral  library*  of  materials  and  then  imaging 
a  target  with  a  known  mix  of  these  materials.  Three  materials  were 
selected  for  the  experiment;  red,  green,  and  yellow  dye.  The  targets 
were  a  known  mix  of  red  and  green  dye.  Yellow  was  produced  where  the 
red  and  green  overlapped  in  the  targets.  The  overlap  is  slight  but 
uncontrolled  and  the  red/green  interaction  is  not  always  in  the  same 
proportion.  It  is  caused  by  saturation  of  the  two  dyes  in  the 
cardboard.  The  spectral  library  was  obtained  by  imaging  separately 
'pure'  red,  green,  and  yellow  dye  on  10x10  inch  cardboards  and 
recording  the  spectra.  These  spectra  constitute  the  library  of 
' spectral  prototypes ' . 

2 . 2  TARGETS 

A  10  x  10  inch  cardboard  was  subdivided  into  64  1J  x  1J  inch 
squares,  8  across  and  8  down.  Each  square  was  alternately  dyed  red 
and  green  so  that  the  cardboard  had  an  even  mix  of  red  and  green  dye, 
32  red  squares  and  32  green  squares.  This  constitutes  a  50/50  percent 
red/green  target.  A  second  target  was  constructed  like  the  first, 
except  that  every  other  green  square  was  dyed  red  so  that  this  target 
had  a  mix  of  48  red  squares  and  16  green  squares.  This  constitutes  a 
75/25  percent  red/green  target.  A  third  target  was  constructed  like 
the  second  one,  except  that  the  red  and  green  dye  are  interchanged  so 
that  this  target  had  16  red  squares  and  48  green  squares.  This 
constitutes  a  25/75  percent  red/green  target.  This  data  set  of 
targets  was  labeled  ’large  square'  targets.  A  second  data  set  was 
constructed  just  like  the  one  above  except  that  the  cardboard  was 
subdivided  into  128  %  x  inch  squares,  16  across  and  16  down. 


This  data  set  was  labeled  'small  square'  targets.  Fig.  1  shows  the 
spectral  library  cardboards  and  the  different  targets. 


2 . 3  SPECTROMETER 

The  instrument  used  for  this  effort  was  a  tripod  mounted/  factory 
calibrated/  Geophysical  Environmental  Research  (GER)  IRIS  MkIV  dual 
field  of  view  spectrometer.  Fig.  2  shows  the  spectrometer  in  a  field 
environment.  Radiometric  values  for  431  discrete  spectral  bands, 
wavelengths  from  350-2500  nanometers  at  5  nanometer  increments,  were 
obtained  for  each  of  the  spectral  prototypes  and  the  targets.  The 
spectrometer's  raw  output  values,  32-bit  hexadecimal  numbers,  were 
downloaded  to  an  IBM-compatible  PC  for  conversion  to  floating  point 
numbers  and  percent  reflectance.  A  Halon  reflectance  standard  was 
used  as  a  reference  for  the  conversion  to  percent  reflectance.  Table 
I  shows  a  sample  of  the  PC  output  data  for  three  of  the  targets.  The 
last  column  of  Table  I  is  a  description  of  the  data  entries.  The  data 
were  collected  in  a  laboratory  with  constant  temperature,  pressure, 
and  lighting  conditions.  It  was  collected  over  two  consecutive  days 
thereby  minimizing  the  probability  of  any  changes  which  could  affect 
the  spectral  information. 

2 . 4  IMAGING  PROCEDURE 

The  spectrometer  was  set  up  on  a  tripod  looking  straight  down,  one 
meter  from  the  floor.  The  checkerboard  target  was  placed  on  the  level 
floor  directly  below  the  spectrometer.  The  board  was  constructed  so 
that  it  would  be  larger  than  the  spectrometer's  circular  field  of 
view.  The  board  was  positioned  so  that  its  geometric  center  was  co¬ 
incident  with  the  spectrometer's  field  of  view  center  so  that  whatever 
portion  of  the  board  was  imaged  by  the  spectrometer,  the  color  dye 
proportion  was  retained.  Fig.  3  shows  how  a  25/75  percent 
checkerboard  was  imaged.  The  field  of  view  is  represented  by 
concentric  circles.  The  portion  of  the  board  inside  each  circle  has 
the  same  red/green  ratio  as  the  whole  board.  This  is  true  as  long  as 
the  targets  are  symmetric  about  the  center,  the  centers  of  the  field 
of  view  and  targets  are  coincident,  and  the  field  of  view  is  within 
the  board  boundaries.  Even  if  the  field  of  view  is  oval  or 
rectangular,  the  red/green  ratio  is  retained  as  long  as  the  above 
three  conditions  hold.  The  '50/50  Red/Green  Large  Squares  Target*  was 
imaged  three  times.  The  second  data  set  is  labeled  'repeated'  and  the 
third  data  set  is  labeled  'repeated  again'.  Other  targets  were 
rotated  45°  for  a  second  image  and  are  labeled  'rotated'. 


FIGUHE1.  SPECTRAL  PHOTOTYPES  AND  TARGETS 


TARGETS:  SMALL  SQUARES 


FIGURE  3.  TARGET  INSTRUMENT  FIELD  OF  VIEW  DIAGRAM 


Checkerboard  =  25/75  Percent  Red. Green  Target 
Concentric  Circles  =  Spectrometer's  Field  of  View 


TABLE  1.  PC  OUTPUT  DATA  SAMPLE 
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File  Name 
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3.  RESULTS 


The  results  obtained  embodied  different  factors,  computer  equipment, 
data,  computations,  and  interpretation. 

3.1  COMPUTER  EQUIPMENT 

The  data  sets  were  processed  on  a  SUN  SPARC  workstation  using  the 
UNIX  operating  system  and  the  X-Windows  windowing  system.  Two 
commercial  software  packages,  Interactive  Data  Language  (IDL)  and 
Mathematica,  were  used  for  the  computations  and  analyses.  'IDL*  is  a 
computing  environment  for  the  interactive  analysis  and  visualization 
of  scientific  and  engineering  data.  It  provides  a  broad  range  of  high 
quality  mathematical  analysis  and  graphical  display  techniques.  The 
graphical  and  statistical  analyses  were  conducted  with  the  aid  of  IDL. 
'Mathematica'  is  a  general  software  system  and  language  intended  for 
matrix  operations,  statistical  analysis,  and  other  mathematical 
applications.  The  matrix  operations  and  mathematical  computations 
were  conducted  with  the  aid  of  Mathematica. 


3.2  DATA  FILTERING 


High  priority  was  given  to  the  data  filtering  process  to  minimize 
the  probability  of  using  false  data.  The  data  suspected  of  being 
contaminated  was  filtered  out.  As  a  consequence,  a  significant  amount 
of  data  were  not  used.  A  plot  of  wavelength  vs.  percent  reflectance 
was  made  of  all  the  data  sets  to  analyze  the  characteristics  of  the 
spectra  and  to  detect  any  anomalies.  Figs.  4-7  show  the  plots  made. 
The  first  few  data  points,  around  350-500  nanometers,  of  all  sets  are 
irregular  spikes  and  seem  to  be  an  idiosyncrasy  of  the  spectrometer 
when  it  first  starts  to  record  data.  These  data  points  are  obvious 
misreadings  since  reflectance  values  greater  than  100  percent  were 
recorded.  These  data  points  were  filtered  from  any  processing.  The 
region  around  800-1100  nanometers  of  all  data  sets  deviates  from  the 
smooth  spectrum  curve  for  no  readily  explanatory  reason  other  than  it 
is  an  'instrument  error'.  Processing  with  or  without  these  data 
points  made  very  little  difference  in  the  computed  solutions.  This  is 
an  indication  that  the  data  were  valid,  however  the  results  which  will 
be  presented  later  do  not  contain  these  data  points.  Data  associated 
with  wavelengths  greater  than  1500  nanometers  seem  to  be  the  same 
regardless  of  the  dye  combination,  probably  because  at  these 
frequencies  the  spectrometer  was  responding  to  the  cardboard  material 
and  not  the  dye.  These  data  points  were  filtered  from  processing. 
Additional  data  points  were  filtered  out  from  individual  data  sets 
which  accounts  for  the  different  number  of  data  points  used  in  the 
various  solutions. 

3 . 3  COMPUTATIONS 

A  computer  routine  was  prepared  to  process  the  collected  data  and 
compute  'coefficients  of  concentrations'  according  to  equation  (4) 
above  and  to  compute  simple  statistics  of  the  results.  Each  data  set 
was  processed  for  a  'two  vector'  solution  where  the  assumption  was 
that  the  mixed  pixel  consisted  of  red  and  green  dye  and  for  a  'three 
vector'  solution  where  the  assumption  was  that  it  consisted  of  red, 
green,  and  yellow  dye. 


These  coefficients  were  converted  to  percent  concentration  (pc)  by. 


pCC^lBi/  £  #*]  100 


Jt-1 


where : 

pc,  =  percent  concentration  for  material  i, 
Si  =  computed  coefficient  for  material  i, 

sk  =  computed  coefficient  for  material  k. 


FIGURE  4 .  PLOTS  OF  SPECTRAL  PROTOTYPES 


N  (_  0X4—— «  0)  t>4->  QCUt) 


X  i_  4X4— <  4)  o+j  <0  C  0  4) 


wave i eng 


25/752  LS  Red/Green  Repeated  PLOT  50/502  LS  Red/Green  PLOT 


x  c.  h  a>  o-p  incut) 


n  u  «t4— <  v  op  ticot) 


500  1000  1500  2000  2500  0  500  1000  1500  2 

wavelength  <nw>  wavelength  <n«> 


500  1000  1500  2000  2500  0  500  1000  1500  2000 

wavelength  <nm)  wavelength  (nm) 


N  (-  4*+--'  4)  OP  QCUO 


*4  (_  4*«— •  4)  OP  <0  C  O  4) 


)  1000  1500  2000  2500  0  500  1000  1500  2000 

wavelength  (nit)  wavelength  (nn) 


Tables  2-7  show  the  results  of  the  various  tests.  The  tables  show  the 
percent  red,  green,  and  yellow,  the  statistics  variance,  standard 
deviation  (std),  mean,  maximum  and  minimum  residuals  for  the  various 
targets.  The  percent  concentration  is  the  most  fundamental  test  of  the 
computed  information  because  if  this  value  is  not  close  to  the  true 
concentration,  the  rest  of  the  statistical  information  is  meaningless. 
As  can  be  seen  from  the  tables,  all  the  computed  concentrations  are 
very  close  to  the  known  concentrations.  It  was  pointed  out  earlier 
that  yellow  was  produced  in  the  targets  where  red  and  green  came 
together.  It  is  not  possible  to  accurately  determine  the  true 
concentration  of  yellow  in  the  targets,  but  all  the  red/green  targets 
have  a  little  yellow.  The  small  square  targets  should  contain  twice 
as  much  yellow  as  the  corresponding  large  ones,  but  only  if  the 
red/green  dyes  overlap  exactly  in  the  same  constant  manner  in  the 
different  targets.  In  the  computations,  the  yellow  dye  in  the  targets 
varied  from  1.56  to  5.16.  Table  8  list-ranks  the  computed  yellow  in 
each  of  the  targets.  The  higher  values,  except  for  one,  belong  to  the 
'small  square'  target  group.  This  seems  to  confirm  the  fact  that 
'small  square*  targets  contain  more  yellow  than  'large  square',  but 
not  in  the  expected  proportion.  The  targets  and  spectral  prototypes 
were  computer  generated  and  printed  with  a  Kodak  XL7700  color  printer. 
Although  the  quality  of  the  printer  is  excellent,  small  variations  in 
paper  quality  and  idiosyncrasies  of  the  printer  can  account  for  some 
of  the  yellow  differences  in  the  targets.  It  was  noted  that  hard  copy 
prints  of  the  generated  spectral  prototypes  contained  slight,  but 
noticeable,  color  variations.  Similar  variation  in  the  targets  could 
exist,  but  in  a  much  smaller  scale  which  would  be  very  hard  to  detect 
and  correct.  Also,  yellow  is  defined  as  a  homogenous  mixture  of  equal 
amounts  of  pure  red  and  pure  green.  The  yellow  spectral  prototype  was 
generated  with  this  definition.  In  the  targets  equal  amounts  of  red 
and  green  produced  yellow,  but  any  deviation  in  the  red/green  ratio 
would  produce  a  shade  of  yellow  whose  spectra  is  different  from  pure 
yellow.  The  model  then  would  not  detect  all  of  these  shades  as 
yellow.  Residual  vectors  were  computed  for  each  data  set  by; 

residual  vector  =  measured  vector  -  computed  vector 

The  computed  vector  is  the  matrix  product  of  the  computed  coefficients 
and  the  spectral  library.  A  large  residual  variance  or  large 
residuals  would  indicate  that  the  data  does  not  fit  the  model,  or  in 
other  words,  that  the  model  is  not  valid.  As  seen  from  the  tables, 
the  residual  variances  vary  from  0.229514  to  2.07936  percent.  These 
low  values  are  indications  of  a  valid  model.  The  size  and 
distribution  of  the  residuals  are  shown  in  Fig.  8-13.  The 
distribution  is  balance,  but  residuals  at  the  higher  frequencies  seem 
to  have  a  different  pattern  than  those  at  the  lower  ones. 


TABLE  2.  PERCENT  CONCENTRATION  RESULTS 


Large  Squares  50/50  Red/Green  Mix 
62  Data  Points  Used 


Red 

Green 

Residual  Variance 
Residual  Standard 
Residual  Mean 
Residual  Max  Value 
Residual  Min  Value 


48.3615% 

51.6385% 

0.399593 

0.632134 

0.117161 

1.328410 

-1.40322 


Large  Squares  50/50  Red/Green  Mix 
62  Data  Points  Used 


Red 

Green 

Yellow 

Residual  Variance 
Residual  Standard 
Residual  Mean 
Residual  Max  Value 
Residual  Min  Value 


47.6452% 

50.7926% 

1.56212% 

0.229514 

0.479076 

-0.0542426 

1.07 

-1.30296 


Large  Squares  25/75  Red/Green  Mix 
71  Data  Points  Used 


Red 

Green 

Residual  Variance 
Residual  Standard 
Residual  Mean 
Residual  Max  Value 
Residual  Min  Value 


27.9018% 

72.0982% 

0.717481 

0.847043 

0.412358 

2.098570 

-2.58418 


Large  Squares  25/75  Red/Green  Mix 
71  Data  Points  Used 


Red 

Green 

Yellow 

Residual  Variance 
Residual  Standard 
Residual  Mean 
Residual  Max  Value 
Residual  Min  Value 


24.9470% 

72.3506% 

2.70239% 

0.349629 

0.591294 

0.042242 

1.204380 

-2.17339 


TABLE  3.  PERCENT  CONCENTRATION  RESULTS 


Large  Squares  75/25  Red/Green  Mix 
73  Data  Points  Used 


Red 

Green 

Residual  Variance 
Residual  Standard 
Residual  Mean 
Residual  Max  Value 
Residual  Min  Value 


75.3445% 

24.6555% 

0.546755 

0.739429 

0.326348 

1.585870 

-2.00696 


Large  Squares  75/25  Red/Green  Mix 
73  Data  Points  Used 


-  73.1665% 

=  24.5850% 

=  2.24847% 

=  0.271439 

-  0.520998 
=  0.0251939 
=  1.576690 
*  -1.59279 

Large  Squares  75/25  Red/Green  Mix  (Rotated) 
72  Data  Points  Used 


Red 

Green 

Yellow 

Residual  Variance 
Residual  Standard 
Residual  Mean 
Residual  Max  Value 
Residual  Min  Value 


Red 

Green 

Residual  Variance 
Residual  Standard 
Residual  Mean 
Residual  Max  Value 
Residual  Min  Value 


75.4534% 

24.5466% 

0.614468 

0.783880 

0.325845 

1.537690 

-1.51769 


Large  Squares  75/25  Red/Green  Mix  (Rotated) 
73  Data  Points  Used 


Red 

Green 

Yellow 

Residual  Variance 
Residual  Standard 
Residual  Mean 
Residual  Max  Value 
Residual  Min  Value 


73.1881% 

24.2567% 

2.55518% 

0.249141 

0.499140 

-0.00771091 

1.45873 

-1.16877 


TABLE  4.  PERCENT  CONCENTRATION  RESULTS 


Large  Squares  50/50  Red/Green  Mix  (Repeated) 
84  Data  Points  Used 


Red 

Green  = 
Residual  Variance  = 
Residual  Standard  = 
Residual  Mean  = 
Residual  Max  Value  = 
Residual  Min  Value  = 


51.1206% 

48.8794% 

0.769353 

0.877128 

0.241784 

2.036880 

-2.04926 


Large  Squares  50/50  Red/Green  Mix  (Repeated) 
70  Data  Points  Used 


Red 

Green 

Yellow 

Residual  Variance 
Residual  Standard 
Residual  Mean 
Residual  Max  Value 
Residual  Min  Value 


48.5735% 

49.0389% 

2.38761% 

0.442796 

0.665429 

0.0265032 

1.792760 

-1.73661 


Large  Squares  50/50  Red/Green  Mix  (Repeated  Again) 
70  Data  Points  Used 


Red 

Green 

Residual  Variance 
Residual  Standard 
Residual  Mean 


48.5496% 

51.4504% 

1.396780 

1.181850 

0.309275 


Residual  Max  Value  =  2.828240 


Residual  Min  Value  =  -2.46116 


Large  Squares  50/50  Red/Green  Mix  (Repeated  Again) 
71  Data  Points  Used 


Red 

Green 

Yellow 

Residual  Variance 
Residual  Standard 
Residual  Mean 
Residual  Max  Value 
Residual  Min  Value 


47.2529% 

48.6534% 

4.09364% 

0.821457 

0.906343 

0.0171588 

3.294860 

-2.66248 


TABLE  5.  PERCENT  CONCENTRATION  RESULTS 


Small  Squares  50/50  Red/Green  Mix 
64  Data  Points  Used 


Red 

Green 

Residual  Variance 
Residual  Standard 
Residual  Mean 
Residual  Max  Value 
Residual  Min  Value 


47.728% 

52.272% 

2.07936 

1.442 

0.488669 

3.32402 

-2.92235 


Small  Squares  50/50  Red/Green  Mix 
65  Data  Points  Used 


Red 

Green 

Yellow 

Residual  Variance 
Residual  Standard 
Residual  Mean 
Residual  Max  Value 
Residual  Min  Value 


45.3574% 

49.5448% 

5.09777% 

0.697811 

0.835351 

-0.0189529 

2.80651 

-1.70586 


Small  Squares  50/50  Red/Green  Mix  (Rotated) 
63  Data  Points  Used 


Red 

Green 

Residual  Variance 
Residual  standard 
Residual  Mean 
Residual  Max  Value 
Residual  Min  Value 


46.4724% 

53.5276% 

1.60674 

1.26757 

0.295009 

3.2525 

-2.16634 


Small  Squares  50/50  Red/Green  Mix  (Rotated) 
60  Data  Points  Used 


Red 

Green 

Yellow 

Residual  Variance 
Residual  Standard 
Residual  Mean 
Residual  Max  Value 
Residual  Min  Value 


45.399% 

50.2944% 

4.30659% 

0.749072 

0.86549 

0.0065667 

2.94102 

-1.56352 


TABLE  6.  PERCENT  CONCENTRATION  RESULTS 


Small  Squares  25/75  Red/Green  Mix 
63  Data  Points  Used 


Red 

Green 

Residual  Variance 
Residual  Standard 
Residual  Mean 
Residual  Max  Value 
Residual  Min  Value 


=  24.1692% 
=  75.8308% 
=  1.95059 
=  1.39664 
=  0.608475 
=  3.1749 
=  -4.60226 


Small  Squares  25/75  Red/Green  Mix 
63  Data  Points  Used 


Red 

Green 

Yellow 

Residual  Variance 
Residual  Standard 
Residual  Mean 
Residual  Max  Value 
Residual  Min  Value 


19.8927% 

74.9423% 

5.16505% 

0.80603 

0.897791 

0.0905302 

1.93587 

-2.78372 


Small  Squares  25/75  Red/Green  Mix 
63  Data  Points  Used 


Red 

Green 

Residual  Variance 
Residual  Standard 
Residual  Mean 
Residual  Max  Value 
Residual  Min  Value 


24.724% 

75.276% 

1.93654 

1.3916 

0.523531 

3.11238 

-3.996 


Small  Squares  25/75  Red/Green  Mix 
63  Data  Points  Used 


Red 

Green 

Yellow 

Residual  Variance 
Residual  Standard 
Residual  Mean 
Residual  Max  Value 
Residual  Min  Value 


21.0989% 

74.4886% 

4.41252% 

1.08177 

1.04008 

0.0763604 

2.05594 

-2.42596 


(Rotated) 


(Rotated) 


TABLE  7.  PERCENT  CONCENTRATION  RESULTS 


Small  Squares  75/25  Red/Green  Mix 
63  Data  Points  Used 


Red 

Green 

Residual  Variance 
Residual  Standard 
Residual  Mean 
Residual  Max  Value 
Residual  Min  Value 


75.964% 

24.036% 

1.90755 

1.38114 

0.399542 

3.64191 

-3.37843 


Small  Squares  75/25  Red/Green  Mix 
63  Data  Points  Used 


Red 

Green 

Yellow 

Residual  Variance 
Residual  Standard 
Residual  Mean 
Residual  Max  Value 
Residual  Min  Value 


75.126% 

21.1961% 

3.67788% 

1.31538 

1.1469 

0.0340338 

3.92794 

-2.0951 


Small  Squares  75/25  Red/Green  Mix  (Rotated) 
63  Data  Points  Used 


Red 

Green 

Residual  Variance 
Residual  Standard 
Residual  Mean 
Residual  Max  Value 
Residual  Min  Value 


77.5492% 

22.4508% 

1.98097 

1.40747 

0.450635 

3.61223 

-3.88359 


Small  Squares  75/25  Red/Green  Mix  (Rotated) 
63  Data  Points  Used 


Red 

Green 

Yellow 

Residual  Variance 
Residual  Standard 
Residual  Mean 
Residual  Max  Value 
Residual  Min  Value 


76.6859% 

19.1925% 

4.12161% 

1.28975 

1.13567 

0.0522504 

3.92398 

-2.48484 


FIGURE  9.  RESIDUAL  PLOTS 

75/25  R/G  LS  Residuals  <rg)  75/25  R/G  LS  Residuols 


FIGURE  11.  RESIDUAL  PLOTS 

50/50  R/G  SS  Residuals  (rg)  50/50  R/G  SS  Residuals  (r 


FIGURE  13.  RESIDUAL  PLOTS 

75/25  R/G  SS  Residuals  (rg)  75/25  R/G  SS  Residuols 
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TABLE  8.  PERCENT  YELLOW  IN  TARGETS 


Rank 

Percent  Yellow 

Red/Green  Combination 

Size  Square 

1 

5.15 

25/75 

Small 

2 

5.09 

50/50 

Small 

3 

4.09 

50/50  (Repeated 
again) 

Large 

4 

4.41 

25/75  (Rotated) 

Small 

5 

4.30 

50/50  (Rotated) 

Small 

6 

4.12 

75/25  (Rotated) 

Small 

7 

3.68 

75/25 

Small 

8 

2.56 

75/25 

Large 

9 

2.70 

25/75 

Large 

10 

2.39 

50/50  (Repeated) 

Large 

11 

2.25 

75/25  (Rotated) 

Large 

12 

1.56 

50/50 

Large 

No  effort  was  made  to  determine  the  various  errors  which  may  influence 
the  results.  Equation  (4)  above  has  an  inherent  error  term,  but  was 
assumed  to  be  zero  for  this  effort.  A  misalignment,  uneven  dyes, 
uneven  checkerboards  and  a  number  of  other  factors  would  cause  errors, 
however  in  light  of  the  analysis  of  the  results  and  statistics,  it 
appears  that  an  error  analysis  at  this  point  would  not  provide  much 
information. 


4.  CONCLUSIONS 

The  results  obtained  show  that  a  linear  estimation  model  and  a 
library  of  spectral  prototypes  are  viable  methods  for  de-mixing 
hyperspectral  mixed  pixel.  The  computed  dye  concentrations  agree  with 
the  known  concentrations  to  a  small  percent  difference.  The  variance 
and  standard  deviation  of  the  residuals  are  indicative  of  a  valid 
model.  The  residuals  themselves  show  no  spectral  bias.  Fig.  8-13 
show  a  plot  of  the  residuals  for  the  various  tests.  In  the  figures 
'LS'  means  large  squares,  'SS'  small  squares,  ’  rg'  red/green  vector, 
and  'rgy'  red /green /yellow  vector. 
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